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ASTIA  AVAIIABIUTY  NOTICt 

QUAllFliS  BiOMSTORt  MAY  OBTAIN  COFIN 
OF  THIS  IBFOBT  FROM  ABTI A. 


ABSTRACT 


During  the  second  year  of  research  under  the  above  contract  the 
analysis  of  the  problem  concerning  the  stability  of  crystal  structures  for 
the  heavy  rare  gas  atoms  has  been  completed.  There  exists  firm  experi¬ 
mental  evidence  that  neon,  argon,  krypton  and  xenon  crystallize  in  the 
face-centered  cubic  (fee)  structure,  whereas  helium,  under  pressure, 
forms  crystals  of  hexagonal  closest  packing  (hep).  According  to  calcul-^ 
ations  based  on  pair- potentials,  however,  the  hexagonal  lattice  should 
be  somewhat  more  stable  for  all  rare  gas  solids  involved. 

This  problem,  which  was  discussed  in  detail  for  the  first  time  in 
the  literature  exactly  twenty  years  ago  (B.  M.  Axilrod  and  E.  Teller, 

J.  Chem.  Phys.  1_1^,  299  (1 943)),  received  much  increased  attention 
since  about  1955.  A  satisfactory  solution,  which  would  at  the  same  time 
constitute  the  first  absolute  stability  proof  for  any  crystalline  structure, 
had,  however,  so  far  not  been  found. 

The  underlying  hypothesis  for  the  analysis  of  this  important  pro¬ 
blem  under  the  present  contract  was  that  the  fcc-structure  becomes  the 
more  stable  one  if  deviations  from  the  principle  of  pair-potentials  bet¬ 
ween  the  atoms  are  taken  into  account.  Since,  in  principle,  there  occur 
many  possible  types  of  such  manybody  interactions,  a  systematic  ana¬ 
lysis  of  these  various  effects  was  undertaken  during  the  first  year  of  the 
contract,  based  on  an  exchange  cluster  expansion  and  making  use  of  a 
perturbation  treatment  (L.  Jansen.  Phys.  Rev.  125,  1798  (1962);  Special 
Technical  and  Scientific  Report  No  2.  July  31,  1961). 

First-order  calculations  were  carried  out  during  the  first  year 
(see  above  references);  the  second  year  was  devoted  to  the  evaluation  of 
second-order  threebody  interactions  between  rare  gas  atoms,  and  to  an 
application  of  the  results  in  first  and  second  orders  to  the  stability  pro¬ 
blem  of  rare  gas  solids. 


To  avoid  difficulties  of  convergence  with  the  use  of  a  multipole  ex¬ 
pansion  for  the  electrostatic  interactions  between  the  atoms,  the  integrals 
were  evaluated  with  the  perturbation  terms  in  unexpanded  form.  Instead, 
for  the  present  type  of  problem  it  was  found  possible  to  use  a  Gaussian- 
effective-electron  model,  in  which  a  single  parameter  characterizes  the 
rare  gas  atom  considered. 

Since  the  first-  and  second-order  threebody  interactions  are  of 
very  short  range,  the  analysis  was  limited  to  triangles  formed  by  a  central 
atom  in  each  lattice  and  two  of  its  twelve  nearest  neighbors.  There  exist 
66  of  such  triangles  for  each  lattice,  57  of  which  are  the  same  between  the 
two  structures,  whereas  9  are  different. 

It  appeared  during  the  analysis  that  the  occurring  many  integrals  of 
sefiond  order  (direct  and  exchange  integrals,  many  of  which  are  of  a  three- 
center  type)  lead  to  long  and  complex  mathematical  problems  for  their 
evaluation.  The  formal  expressions,  being  considerably  more  complicat¬ 
ed  than  those  occurring  in  first  order,  require  high  numerical  precision 
to  render  the  final  result  reliable.  First  essays,  with  electronic  comput¬ 
ation,  proved  much  too  costly  and  not  sufficiently  accurate  for  this  type 
of  problem.  . 

It  was  therefore  necessary  to  look  for  special  methods,  based  on 
integral-transformation  followed  by  asymptotic  series  expansions.  Fortu¬ 
nately,  and  to  our  surprise,  it  appeared  possible  to  find  excellent  asymp¬ 
totic  expansions  for  the  five  basic  integral  types  involved  in  the  problem. 

The  precision  attained  in  this  way  was  much  more  satisfactory  than  that 
given  by  electronic  computation,  especially  for  large  values  of  the  variables. 

In  this  way,  it  was  possible  to  evaluate  the  final  result  with  required 
precision  and  to  express  it  in  analytic  form;  the  expression  is  of  the  type 
of  a  sum  of  exponentially  decreasing  functions  of  the  three  sides  of  the  tri¬ 
angle  considered,  plus  a  sum  of  inverse  powers  of  these  distances. 
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We  considered,  as  in  the  first-order  calculations,  the  nine  different 
triangles  among  the  sixty-six  in  each  lattice,  formed  by  a  central  atom 
and  two  of  its  twelve  nearest  neighbors.  The  results  can  be  expressed  in 
terms  of  the  dimensionless  parameter  p  R.  where  P  is  the  characteristic 
parameter  of  the  Gaussian  charge  distribution,  and  R  the  nearest-neighbor 
distance  in  the  lattice,  and  the  angle  d  between  the  two  sides  of  the  triangle 
which  originate  at  the  central  atom. 

Detailed  results  for  argon  were  obtained  for  0  *  60°,  and  for  the 
critical  120°-  and  180°-configurations  of  the  fcc-structure;  approximate 
values  were  already  determined  for  the  other  important  triangles  0  =110° 
and  0  =  146°  of  the  hcp-lattice.  In  addition,  the  behavior  of  the  relative 
second-order  threebody  component,  aE2/E2^°^,  where  A  Eg  is  the  three- 
body  contribution,  the  pair-potential  second-order  energy  for  the 

atoms  on  the  triangle,  as  a  function  of  8  for  the  other  heavy  rare  gas  atoms 
was  established. 

The  first-order  results,  obtained  previously,  were  then  reexpressed 
in  terms  of  0  and  PR  for  the  same  triangles,  and  the  corresponding  relative 
first-order  threebody  component  aEj/E^^°^  compared  with  its  analogue  in 
second  order.  The  sensational  result  of  this  comparison  is  that  the  curves 
for  aE^/E^^°^  and  AEo/En^°^  are  practically  one  and  the  same  function; 
for  argon  it  starts  at  -20  percent  for  0  =  60°,  increases  very  rapidly  until 
8^  120°  and  is  almost  constant  from  0  *  120°  to  0  »  180°. 

We  have  also  demonstrated  that  this  type  of  curve  is  unambiguously 
in  favor  of  the  hcp-lattice  as  far  as  aEj/Ej^°^  is  concerned.  On  the  other 
hand,  the  curve  stabilizes  the  fcc-lattice  in  second  order,  since  the  de¬ 
nominator  Eg^®^  is  negative.  Since  the  lattice  energy  itself  is  a  negative 
quantity,  the  resulting  curve  favors  the  face-centered  cubic  configuration. 

A  more  accurate  comparison  for  the  types  of  triangles  mentioned  above 
yielded  aE^  ^  -2  aE^,  which  is,  remarkably  enough,  just  the  same  ratio 
as  there  exists  between  Eg^®^  and  Ej/°\ 
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This  effect  has  some  of  the  general  characteristics  of  the  Teller- 
Axilrod  triple -dipole  (third-order)  contribution,  but  it  gives  rise  to  a 
much  larger  difference  between  the  fee-  and  hcp-structures.  The  large 
difference  derives  from  the  very  special  behavior  of  the  or 

A  Eg/ Eg curve  in  the  neighborhood  of  0  =  110°  to  0  »  120°. 

In  addition  to  explaining  the  stability  of  the  face-centered  cubic 
configuration  for  solids  of  the  heavy  rare  gases,  the  theory  predicts,  in 
particular,  that  the  interaction  between  atoms  in  an  equilateral  triangular 
array  will  be  20  percent  less  negative  than  the  sum  of  its  pair-potentials. 
In  this  connection,  a  very  recent  investigation  of  the  specific  heat  of 
solid  argon  (at  constant  volume)  by  Foreman  and  Lidiard  (Phil.  Mag.  ,  in 
the  press)  is  of  great  interest.  For  their  calculation,  these  authors  used 
the  anharmonic  Einstein  model;  from  the  results  they  estimated  the  free 
energy  of  formation  of  a  vacancy.  There  is  a  discrepancy  between  these 
results  and  experimental  values.  Foreman  and  Lidiard  suggest  after  care 
ful  analysis  that  the  discrepancy  can  be  due  to  the  effect  of  threebody  in¬ 
teractions  and  estimate  a  relative  effect  of  about  -25  percent  for  an  equi¬ 
lateral  triangle  of  argon  atoms.  Our  theoretical  result  is  in  agreement 
with  their  estimate. 

It  may  possibly  also  be  anticipated  that  the  unexplained  difference 
in  crystal  structures  occurring  among  the  alkali  halides  (for  example, 
NaCl  and  KCl  on  one  hand,  CsCl  on  the  other  hand)  will  find  a  similar  ex¬ 
planation.  We  note  that  the  corresponding  ions  are  isoelect ronic  with  the 
rare  gas  atoms.  Preliminary  investigations  are  in  progress  on  this  sub¬ 
ject. 


We  have  also  initiated  a  band  theoretical  calculation  for  the  same 
problem,  in  which  the  influence  of  symmetry  on  the  positions  and  multi¬ 
plicities  of  the  s-  and  p-bands  of  solid  argon  in  the  fee-  and  hep- configur¬ 
ations  is  studied.  In  this  analysis,  extensive  use  is  made  of  group-theo¬ 
retical  considerations.  Again,  the  first  result  has  been  that  the  difference 


•  5  - 


in  lattice  energy  between  the  two  configurations  must  be  due  to  the  effect 
of  manybody  interactions.  The  various  threebody  terms,  which  are  now  of 
the  type 


j’u*j(r)  V(r  -  R^)  u^  (r  -  Rp)  dx  . 


where  u^,  denote  3s-  or  3p-functions,  V(r  -  R,^)  the  atomic  potential 
of  an  atom  with  nucleus  centered  at  R^,  and  where  p  ^  V,  and  p  ,vf  0, 
are  being  evaluated. 

This  component  of  the  research  must  be  seen  in  connection  with  the 
general  problem  concerning  the  effect  of  (slight)  changes  in  symmetry  of 
a  solid  system  on  its  band  structure;  the  analysis  is  being  undertaken  in 
conjunction  with  a  different  project.  Results  have  already  been  obtained 
for  a  slightly  distorted  fcc-lattice  (2  atoms  per  unit  cell  instead  of  one); 
it  is  envisaged  that  also  the  influence  of  deviations  from  the  ideal  c/a- 
ratio  in  fcc-lattices,  deviations  which  normally  occur  in  nature,  on  the 
band  structure  of  molecular  solids  will  be  analyzed. 


OUTLINE  OF  RESEARCH  PROJECT 

Concerning  properties  of  materials  in  general,  it  is  nowadays  re¬ 
cognized  that  no  essential  progress  is  possible  without  extensive  funda¬ 
mental  research.  In  the  past,  most  of  the  theoretical  investigations  in  the 
field  of  properties  of  solid  materials  have  been  directed  towards  an  ex¬ 
planation  of  phenomena  associated  with  metallic  substances,  of  semi-conr 
ductors,  and,  to  a  lesser  extent,  of  ionic  crystals.  Compared  with  this 
group  of  materials,  the  theory  of  molecular  solids  (insulators)  has  re¬ 
ceived  but  little  attention,  and  even  the  most  striking  facts  eidiibited  by 


such  materials  often  lack  a  satisfactory  explanation.  It  has  becoAie  more 
and  more  urgent  to  try  to  fill  this  gap.  since  molecular  solids  have  become 
of  increasing  importance  also  from  a  practical  point  of  view.  As  an  example 
it  might  be  mentioned  that  the  stability  of  the  close-packed  cubic  crystal 
structure  of  the  heavy  rare  gases  (neon,  argon,  krypton  and  xenon)  cannot 
be  explained  on  the  basis  of  existing  theories.  One  finds  that  the  hexagonal 
close-packed  crystal  form  (the  helium  structure)  should  be  more  stable  than 
the  cubic  structure,  and  this  result  is  remarkably  insensitive  with  respect 
to  the  precise  form  of  intermolecular  potentials.  Concerning  this  problem 
the  present  investigator  has  suggested  that  an  explanation  may  be  found  on 
the  basis  of  the  effect  of  manybody  forces  of  the  van  der  Waals-type  (second- 
order  of  perturbation  theory)  on  the  stability.  Preliminary  calculations  have 
shown  that  this  explanation  may  be  the  correct  one,  and  that  present  theories 
of  molecular  solids  need  considerable  modification. 

The  crystals  of  nitrogen  and  carbon  monoxide  are  also  cubic  close- 
packed  at  the  lowest  temperatures,  but  show  a  transition  to  the  hexagonal 
closest  packing  at  35°K  and  61°K,  respectively.  The  present  investigator 
has  suggested  that  the  stability  of  these  crystals  and  the  transition  mecha¬ 
nism  may  involve  orientational  interactions  (e.  g.  those  due  to  permanent 
electric  quadrupole  moments)  in  the  solids.  The  cubic  structure  is  favored 
by  quadrupole  interactions  in  the  solids,  whereas  the  hexagonal  structure  is 
preferred  by  anisotropic  components  of  the  van  der  Waals  forces.  As  the 
temperature  increases  and  herewith  the  rotational  amplitudes  of  the  mole¬ 
cular  axes,  the  quadrupole  forces,  being  of  first  order,  die  out  much  more 
quickly  that  the  (second-order)  van  der  Waals  forces.  The  results  obtained 
so  far  suggest  strongly  that  this  explanation  may  be  the  correct  one,  although 
a  quantitative  theory  is  still  lacking. 

Molecular  solids  have  recently  become  of  importance  also  as  materials 
in  which  free  radicals  can  be  trapped  and  stored.  The  results  obtained  in  this 
field  provide  an  invaluable  experimental  background  for  the  development  of  a 


crystal  field  theory  of  molecular  solids.  This  investigator  has  been  asso¬ 
ciated  with  the  Free  Radicals  Research  Group  at  the  National  Bureau  of 
Standards  in  Washington,  D.  C. ,  and  has  performed  some  preliminary  cal¬ 
culations  which  seem  to  indicate  that  present  theories  of  molecular  solids 
are  in  need  of  substantial  improvement  in  order  to  yield  satisfactory  agree¬ 
ment  with  experimental  data. 


SUMMARY  OF  RESULTS  OBTAINED  DURING  SECOND  YEAR 

The  results  obtained  during  the  second  year  under  the  project  may 
be  divided  into  three  categories: 

I,  Completion  of  the  analysis  concerning  absolute  stability  of  the  cubic 
configurations  of  heavy  rare  gas  solids,  following  the  general  classific¬ 
ation  of  manybody  interactions  developed  during  the  first  year  (Phys.  Rev. 

125,  1798  (1962)).  The  theory  cannot  only  explain  this  stability,  it  is  also 
in  agreement  with  recent  experimental  results  concerning  the  free  energy 

of  formation  of  a  vacancy  in  solid  argon  (A.  J.  E.  Foreman  and  A.  B.  Lidiard, 
Phil.  Mag. ,  in  the  press). 

Two  manuscripts  containing  the  results  of  this  investigation  ("Absolute 
Stability  of  Cubic  Crystal  Structures  of  Heavy  Rare  Gas  Atoms"  I  and  II) 
will  be  submitted  for  publication  in  Physics  Letters  before  Februaryr  1.  Their 
texts  are  given  in  original  form  in  the  following  section  ANALYSIS  OF  RE¬ 
SEARCH  PERFORMED  of  this  Report. 

II.  Development  of  asymptotic  series  expansions  for  integrals  of  three- 
center  type.  The  corresponding  methods  and  results  for  various  triangular 
configurations  of  rare  gas  atoms  will  be  published  elsewhere  and  also  pre¬ 
pared  for  the  Special  Technical  and  Scientific  Report  No.  6. 


in.  A  band  theoretical  analysis  of  the  argon  stability  problem.  We  await 
here  the  evaluation  of  three- center  integrals.  Dr  Michael  H.  Boon  is  con¬ 
ducting  this  part  of  the  research,  in  collaboration  with  Dr  Erminio  Lombardi, 
a  guest  from  the  Institute  of  Technology,  Milano,  Italy,  Although  Dr  Lombardi 
is  engaged  in  research  under  the  project^  his  expenditures  are  covered  by  in¬ 
dependent  sources. 

A  Special  Technical  and  Scientific  Report  No. 5  (May  31,  1962)  was 
issued  containing  some  preliminary  results  on  the  evaluation  of  three-center 
integrals. 


-  8  - 


ANALYSIS  OF 


RESEARCH  PERFORMED 


GENERAL 

The  following  sections  contain  a  discussion  of  results  obtained  on 
the  problem  of  stability  of  crystals  of  the  heavy  rare  gas  atoms.  For 
simplicity,  the  texts  of  two  manuscripts  with  the  title  ABSOLUTE  STABI¬ 
LITY  OF  CUBIC  CRYSTAL  STRUCTURES  OF  HEAVY  RARE  GAS  ATOMS, 
I  and  II,  are  given  in  their  original  forms,  whereby  the  page-numbering 
of  these  manuscripts  has  been  retained.  The  texts  have  been  prepared 
for  publication  in  Physics  Letters  (Amsterdam),  since  it  was  felt  that  the 
results  warrant  quick  publication. 

Since  the  band  theoretical  analysis  of  the  stability  problem  for  solid 
argon  has  not  yet  been  completed,  the  results  obtained  so  far  on  Uiis  part 
of  the  research  are  not  included  in  this  Report. 

A  brief  discussion  of  the  results  on  the  stability  problem  has  been 
presented  in  the  first  section  of  the  present  Report;  we  will,  therefore, 
not  abstract  these  results  here. 


ABSOLUTE  STABILITY  OF  CUBIC  CRYSTAL  STRUCTURES  OF 
HEAVY  RARE  GAS  ATOMS  I* 


Laurens  Jansen 

Battelle  Memorial  Institute.  International  Division, 
Geneva,  Switzerland 

(Received  ) 


Twenty  years  ago^^  the  problem  concerning  the  absolute  stability 
of  the  observed  face-centered  cubic  crystal  structures  of  neon,  argon, 
krypton  and  xenon  was  for  the  first  time  discussed  in  detail  in  the 
literature.  In  conjunction  with  the  pioneering  work  by  Born  and  colla¬ 
borators^^  on  the  mechanical  stability  of  crystals,  this  also  marked  the 
first  realistic  attempt  to  provide  an  explanation  for  the  absolute  stabi¬ 
lity  of  any  crystalline  substance. 

The  problem  arose  since  it  was  found  that  all  other  available  experi¬ 
mental  data  on  crystals  of  neon,  argon,  krypton  and  xenon  could  be  explain- 
ed  in  a  satisfactory  manner  on  the  basis  of  pair- potentials  between  the 
atoms.  However,  the  static  lattice  energy,  written  as  a  sum  of  inter¬ 
actions  between  pairs  of  atoms,  favors  the  hexagonal  close-packed  (hep) 
structure  instead  of  the  experimentally  found  face-centered  cubic  (fee) 
configuration.  The  difference  in  energy  between  these  two  lattices  is  very 
small,  only  of  the  order  of  one  hundredth  of  one  percent,  but  is  remark¬ 
ably  stable  with  respect  to  admissible  changes  in  the  potential  between 
two  atoms 

This  research  has  been  made  possible  through  the  support  and  sponsor¬ 
ship  of  the  U.  S.  Department  of  Army,  through  its  European  Research 
Office. 
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Careful  considerations  of  the  influence  of  zero- point  energy^^®^, 
and  of  the  possibility  for  thermal  transitions  between  the  two  structures 
by  Barron  and  Domb^^  yielded  the  same  negative  results  :  the  hcp-stru- 

4 

ture  is  somewhat  more  stable  than  the  fcc-configuration.  Only  He  ,  crys¬ 
tallizing  under  pressure  in  the  hcp-structure,  appeared  to  behave  normal¬ 
ly. 

*  There  are  two  categories  of  possible  solutions  for  this  problem  : 

a)  the  interactions  between  rare  gas  atoms  contain  non- negligible 
manybody  components,  i.  e.  the  crystal  energy  cannot  be  written 
as  an  additive  sum -over-pairs. 

b)  the  interaction  between  two  atoms  is  angular  dependent. 

Practically  all  proposed  solutions  have,  explicitly  or  implicitly, 
been  based  on  assumption  a).  Only  Cuthbert  and  Linnett^^  suggested  an 
explanation  of  type  b)  ;  according  to  calculations  by  Linnett  and  P6e^*^^ 
the  most  probable  configuration  of  the  group  of  eight  electrons  in  heavy 
rare  gas  atoms  is  as  four  pairs  at  the  corners  of  a  regular  tetrahedron, 
the  nucleus  being  at  the  center.  This  leads  to  an  angular  dependent  electro¬ 
static  interaction  between  two  such  tetrahedra,  which  favors  the  cubic  lat'. 
tice.  No  quantitative  results  were  obtained. 

On  the  other  hand,  the  presumedly  very  high  stability  of  the  closed- 
shell  configurations  of  rare  gas  atoms  led  to  the  supposition  that  many¬ 
body  interactions  between  such  atoms  could  only  be  a  higher-order  per¬ 
turbation  effect  (starting  from  third  order),  or  that  they  could  be  im¬ 
portant  only  at  very  small  interatomic  distances,  where  first-order  inter¬ 
actions  prevail^^^.  The  first  attempt,  by  Axilrod^^\  was  a  third-order 
perturbation  calculation  for  three  non-overlapping  rare  gas  atoms,  i.  e. 
a  straightforward  extension  of  London's  derivation  for  second-order  di- 
pole-dipole  interactions.  This  so-called  triple-dipole  effect  decreases  the 


* 
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attraction,  compared  with  an  additive  sum-over-pairs,  for  an  equilateral 
configuration  of  atoms,  and  Increases  the  attraction  for  a  linear  array.  The 
net  effect,  summed  over  the  hep-  and  fcc-lattlces,  favors  the  cubic  con- 
fi^ration,  but  Its  magnitude  Is  too  small  to  overcome  the  small  palr-potent- 
lal  difference  between  the  two  structures. 

First-order  calculations  for  three  helium  atoms  were  carried  out  by 
P,  Rosen^^^  and  Shostak^^^  Although  helium  is  of  no  direct  interest  as  far 
as  its  crystal  structure  is  concerned,  the  sign  of  this  first-order  threebody 
effect  may  be  noted  :  repulsion  between  three  atoms  is  weaker  for  an  equi¬ 
lateral  triangle,  and  stronger  for  a  linear  array,  than  that  calculated  on  the 
basis  of  an  additive  sum-over-pairs,  i.  e.  the  relative  first-order  effect  has 
the  same  sign  as  the  third-order  triple-dipole  contribution. 

A  variety  of  different- type  manybody  interactions  has  further  been  re¬ 
ported  in  the  literature,  based  partly  on  electrostatic  effects,  in  which  over¬ 
lap  of  charge  clouds  is  treated  classically^^^®^^®^  or  based  on  the  Drude- 

171-211  ♦ 

model  of  harmonic  oscillators  for  the  atoms  '  '  with  dipole- interactions  , 

Finally,  it  was  shown  by  one  of  us  and  R.  T.  McGinnies^^^^^^  that  also 
second-order  interactions  (i.  e.  of  the  usual  van  der  Waals-type)  are  not  ad¬ 
ditive  if  exchange  effects  are  taken  into  account.  The  calculations  were 
carried  out  for  dipole-dipole  and  dipole-quadrupole  components;  a  tendency 
towards  stabilisation  of  the  face-centered  cubic  configuration  was  observed. 

All  effects  discussed  above  suffer  from  one,  or  several,  of  three  main 
defects  which  make  them  unsuitable  as  a  key  to  the  explanation  of  the  ob¬ 
served  stability  of  the  face- centered  cubic  structures  for  neon,  argon,  kryp¬ 
ton  and  xenon  : 


We  mention  also  that  R.  S.  Mulliken  has  suggested  that  polarization  forces 
of  a  "charge-transfer"  type  may  have  an  effect  on  the  stability  of  rare  gas 
crystals.  (R.  S.  Mulliken,  "The  Interaction  of  Electron  Donors  and  Accep¬ 
tors",  Confe^nce  on  "Quantum  Mechanical  Methods  in  Valence  Theory", 
published  by  Oi'f ice  of  Navsd  Research,  Washington,  D.  C.  i^ept.  1951,  p. 

1 08. )  To  this  suthoir's  knowledge  the  above  suggestion  has  never  been  work¬ 
ed  out. 
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1)  they  concern  systems  of  little  interest  (helium); 

2)  overlap  between  neighboring  charge  clouds  is  either 
neglected  or  treated  classically; 

3)  a  multipole  expansion  is  used  for  the  electrostatic  inter¬ 
actions. 

At  small  distances  between  atoms  (those  between  nearest  neighbors 
in  the  crystal)  the  interactions  are  mainly  of  exchange-type.  Further,  the 
use  of  a  multipole  expansion  for  the  electrostatic  interactions  becomes  in¬ 
creasingly  less  accurate  as  the  distances  between  the  atoms  decrease.  In 
fact,  such  series  usually  do  noi  converge  at  all^^^  Exchange  effects 
normally  reverse  the  sign  of  contributions  based  on  classical  models; 
since  the  sign  of  any  effect  is  of  decisive  importance  for  stability  consider¬ 
ations,  critical  errors  maybe  introduced  in  electrostatic  models. 

To  find  a  satisfactory  solution  for  the  stability  problem  a  model 
should  be  designed  which  avoids  the  defects  mentioned  earlier,  which  is 
sufficiently  simple  for  numerical  calculations,  which,  on  the  other  hand, 
retains  the  possibly  important  factors  characterizing  the  problem.  Such  a 
model  may  be  devised  as  follows.  We  note  that  stability  of  the  fcc-struc- 
tures  is  a  phenomenon  common  to  solid  neon,  argon,  krypton  and  xenon. 
Therefore,  the  explanation  can  depend  only  on  some  general  parameter 
characteristic  for  the  wave  functions  of  these  atoms,  not  on  the  detailed 
analytic  form  of  the  wave  functions.  We  have  here  the  choice  between  de¬ 
viation  from  spherical  symmetry  and  the  dimension  of  an  atom,  and  we 
choose  the  latter  parameter.  With  this  choice,  the  hcp-structure  is  sup¬ 
posed  to  be  the  more  stable  one  for  point-atoms,  and  for  those  solids  of 
spherical  atoms  for  which  the  ratio  between  atomic  dimension  and  nearest- 
neighbor  distance  is  smaller  than  a  certain  critical  value,  to  be  deter¬ 
mined.  For  heavier  atoms  the  fcc-structure  should  become  more  and 
more  stable. 


Further  simplifications  maybe  made.  Explicit  second-order  contri¬ 
butions  of  threebody  interactions  between  (helium  or)  neon  atoms  showed^ 
that  multiple  exchange  of  electrons  between  the  same  two  atoms  is  not  im¬ 
portant  for  solids  of  these  atoms.  In  addition,  coupling  of  inter-  and  intra- 
atomic  exchange  gives  only  a  very  small  contribution  to  the  threebody  ex¬ 
change  energy.  Since  neon  crystallizes  already  in  the  fcc-configuration, 
we  conclude  that  also  for  the  heavier  atoms  multiple  exchange  (between 
the  same  two  atoms)  and  coupled  exchange  can  be  neglected. 

We  are  now  left  with  a  sum  (average)  of  single- exchange  interactions 
between  different  atoms.  This  sum  is  replaced  by  an  exchange  interaction 
between  effective  electrons,  one  on  each  atom.  The  problem  becomes  thus 
formally  the  same  as  that  concerning  hydrogen  atoms  with  parallel  spins. 

A  convenient  choice  for  the  charge  distribution  of  an  effective  electron  is  a 
Gaussian  function 


P  (r)  -  exp  (-  p^r^), 

where  r  is  the  distance  from  the  effective  electron  to  its  nucleus  and  where 
p  is  a  parameter  which  measures  the  width  of  the  charge  distribution.  It 
is  possible  to  determine  the  values  of  p  empirically  by  calculating  van  der 
Waals  interactions  between  two  such  Gaussian  atoms  and  equating  the  re¬ 
sult  to  the  attractive  part  of  an  empirical  potential  function^^\ 

The  threebody  effect  will  then  depend  only  on  the  configuration  of  a 

specific  triangle  considered  and  on  the  value  of  the  dimensionless  parameter 

P  R,  where  R  is  the  nearest-neighbor  distance  between  the  atoms.  For 

solid  neon,  argon,  krypton  and  xenon  we  find  PR  »3.44,  2. 10  and  2.  00, 

2 

respectively.  The  energy  of  interaction  is  expressed  in  units  pe  ;  for  argon, 

2 

for  example,  we  find  that  p  e  «  lOeV,  i,  e. ,  somewhat  smaller  than  its  first 
ionization  potential  (  P  «  0.  623  R  «  3.  84  %.).  For  solid  helium  we  have 
P  R  B  5. 75,  but  validity  of  a  Gaussian  model  for  this  atom  is  very  doubtful. 


-  6  - 


We  should  therefore  find  that  above  a  certain  critical  value  of  pR  the 
hcp-structure  is  favored,  whereas  below  this  value  the  fcc-lattice  becomes 
the  more  stable  configuration.  To  carry  out  the  calculations,  we  note  that 
the  sublimation  energy  of  rare  gas  crystals  is  of  the  order  of  one  hundred 
times  smaller  than  those  of  ionic  crystals  or  of  metals.  A  perturbation  treat¬ 
ment  is,  therefore,  appropriate,  whereby  it  is  supposed  that  only  first-  and 
second-order  effects  must  be  included. 

It  is  further  supposed,  since  exchange  effects  decrease  very  rapidly 
with  interatomic  distances,  that  the  main  effect  which  stabilizes  the  fcc- 
structure  lies  already  in  a  study  of  the  difference  between  triplet-inter- 
actions  involving  a  central  atom  and  any  two  of  its  twelve  nearest  neighbors 
in  the  two  configurations.  As  far  as  pair- interactions  are  concerned  the  two 
structures  are  identical  up  to  the  third  shell  of  neighbors.  However,  if  tri¬ 
angular  configurations  are  taken  into  account,  then  we  count  66  possible  tri¬ 
angles  in  each  structure,  formed  by  a  central  atom  and  two  of  its  nearest 
neighbors.  Of  these  66  triangles,  57  are  the  same  between  the  two  lattices, 
but  9  are  different.  These  9  different  triangles  are  listed  in  Table  I  below; 
a,  b,  and  c  denote  the  three  sides  of  the  triangle  in  units  of  nearest-neighbor 
distance.  We  also  list  values  for  the  angle  9  between  the  two  sides  originating 
at  the  central  atom. 


Table  I ;  The  nine  different  triangles  a,  b,  c  with  a  »  b  =  1,  between  the 
hep-  and  fcc-structures. 


2  2  2 

2  .^2  2 

a  ,  b  ,  c 

9 

No 

a  ,  b  ,  c 

9 

No 

hep 

• 

00 

110° 

3 

1,  1.  11/3 

146° 

6 

fee 

CO 

% 

120° 

6 

1.  1.  4 

180° 

3 
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These  different  triangles  are  all  formed  by  a  central  atom,  one 
atom  on  a  triangle  above  the  horizontal  hexagonal  layer,  and  one  atom  on 
a  triangle  below  this  layer.  Corresponding  atoms  on  the  triangles  above 
and  below  are  just  above  each  other  in  the  hcp-lattice;  in  the  fee- configur¬ 
ation  one  of  the  triangles  is  rotated  through  60*^  about  the  axis  through  its 
center  and  perpendicular  to  the  horizontal  plane. 

A  perturbation  calculation  is  now  carried  out  in  first  and  second 
orders  for  the  interaction  energy  between  three  Gaussian  atoms  in  a  spe¬ 
cific  triangular  configuration.  Different  values  of  the  parameter  p  repre¬ 
sent  the  different  rare  gas  atoms.  The  use  of  a  multipole  expansion  is 
avoided  by  retaining  the  electrostatic  interactions  in  their  exact  form. 

First-order  calculation.  The  correction  to  first-order  (repulsive) 
interactions  between  three  atoms  in  different  triangular  configurations 
has  been  published  recently^^\  These  values  have  been  supplemented  by 
calculating  also  the  correction  for  the  critical  110°-  and  146°-triangle8 
of  the  hcp-lattice.  The  result  is  (cf.  Fig.  1  of  ref.  27)  that  the  146°- 
curve,  as  a  function  of  PR,  lies  about  half-way  between  the  120°-  and 
180°-curves,  and  that  the  effect  for  a  110°-configuration  is  very  nearly 
zero,  independent  of  the  value  for  PR.  In  the  following  Fig.  1  we  plot 
the  relative  threebody  component  of  the  first-order  energy,  AE^/Ej, 
as  a  function  of  the  angle  6,  for  argon  (PR  «  2.  4)  and  xenon  (PR  >  2.  0). 
Note  that  the  sign  of  this  effect  is  the  same  as  that  found  for  three  helium 
atoms  by  P.  Rosen^^^  and  Shostak^^^ 


(Fig.  1) 


Fig.  1 

Relative  first-order,  threebody  interaction 
for  triangles  a  *  b  =1,  c  (units  of  nearest-neighbor 
distance)  of  argon  and  xenon  atoms,  as  a  function  of  the 
angle  6  between  the  sides  a  and  b.  The  symbols  h  and  ^ 
c  denote  the  hep-  and  fcc-lattice,  respectively;  the  num¬ 
ber  of  triangles  at  a  specific  angle  6  is  given  in  parantheses. 
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The  curve  for  as  in  between  those  for  argon  • 

and  xenon,  whereas  the  rt.  der  effect  for  neon  shows  very 

nearly  the  same  behavior  as  ^e  other  rare  gas  atoms  for  6  <  110°; 

for  larger  values  of  0  it  stays  muv  .i  closer  to  the  horizontal  axis. 

It  is  easy  to  show  that  a  typical  curve  as  given  in  Fig.  1  stabilizes 
the  hexagonal  close-packed  configuration.  To  see  this  we  have  indicated 
in  the  figure  the  important  values  0  =  110°  and  146°  (hep)  and  0  =  120° 
and  180°  (fee),  supplied  the  corresponding  points  of  the  curve  with  in¬ 
dices  h  (hexagonal)  and  c  (cubic)  and  written  in  parentheses  the  number 
of  triangles  occurring  for  that  value  of  0  in  the  corresponding  structure. 

To  make  an  accurate  comparison  between  the  two  lattices,  we 
have  to  compare  the  values  of  AE^  for  a  fixed  value  of  itself,  for 
example,  at  120°.  This  transformation  will  flatten  the  curve  between 
120°  and  180°,  leave  the  zero-value  at  110°  unchanged  and  lower  the 
value  somewhat  for  an  equilateral  triangle.  This,  however,  does  not 
change  the  argument  in  any  appreciable  way. 

To  a  very  good  approximation,  the  curve  for  aE^/E^  (0  =  120°) 
increases  (slowly  and)  linearly  with  0  between  120°  and  180°.  Let  us 
indicate  its  value  at  146°  by  X,  at  180°  by  X  +  a .  Then  at  120°  the  value 
is  very  nearly  X  -  a  .  At  1 1 0°  we  have  a  contribution  which  we  denote  by  Y, 
with  Y  »  0.  The  comparison  yields  ; 

• 

[Ej(fcc)  -  Ej(hcp)]/Ej-  6pc  -  a)  +  3pc  +a)  -  (6X  +  3Y)  -  3[(X  -  a)  -  Y]  >0. 

Since  E^  is  positive,  this  means  that  the  fcc-lattice  has  a  higher  (less  ne?- 
gative)  energy  than  the  hep- configuration,  so  that  the  hexagonal  close- 
packed  lattice  is  stabilized  by  first-order  threebody  interactions.  This 
difference  is  an  order  of  magnitude  larger  than  that  in  static  lattice  ener¬ 
gy  based  on  pair-potentials;  it  has,  therefore,  the  desired  magnitude,  but 
the  wrong  sign. 
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The  peculiar  behaviour  of  the  AEj/E^-curve  (Fig.  1)  is  to  be  borne 
in  mind;  it  will  be  of  particular  importance  for  the  second-order  calculation 
which  we  discuss  in  the  following  communication.  The  first-  and  second- 
order  effects  will  then  be  combined  to  determine  the  total  effect  of  threebody 
interactions,  in  the  first  two  orders  of  perturbation  theory,  on  the  stability 
of  the  fee-configurations  for  solids  of  the  heavy  rare  gas  atoms. 
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There  is  a  good  experimental  evidence  that  neon,  argon,  krypton 
and  xenon  crystallize  in  a  face-centered  cubic  (fcc)-lattice,  whereas 
helium,  under  pressure,  forms  crystals  of  hexagonal  closest  packing 
(hep).  This  preference  of  heavy  rare  gas  solids  for  the  cubic  configur¬ 
ation  has  been  the  subject  of  many  investigations,  but  a  satisfactory 
solution  has  so  far  not  been  found. 

In  a  previous  communication^^  this  problem  has  been  outlined. 

It  was  supposed  that  short-range  manybody  interactions  in  these  solids 
are  significant  and  that  the  difference  between  these  effects  in  the  two 
close- packed  structures  (fee  and  hep)  may  be  sufficiently  large  to  over¬ 
come  the  small  difference  in  static  lattice  energy  based  on  pair-potent¬ 
ials.  which  favors  the  hcp-lattice. 


This  research  has  been  ma.de  possible  through  the  support  and  spon¬ 
sorship  of  the  U.  S.  Department  of  Army,  through  its  European  Re¬ 
search  Office. 
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First-order  perturbation  results,  obtained  on  the  basis  of  a  Gaussian 
effective-electron  model,  for  threebody  interactions  were  seen  to  be  unam¬ 
biguously  in  favor  of  the  hep- structure.  The  triangles  considered  are  formed 
by  a  central  atom  and  any  two  of  its  twelve  nearest  neighbors  in  the  two  lat¬ 
tices.  The  remaining  task  is  to  calculate  the  second-order  threebody  inter¬ 
actions  between  heavy  rare  gas  atoms,  and  then  to  combine  the  first-  and 
second-order  effects  to  provide  a  definite  answer  to  the  stability  problem. 
Since  rare  gas  crystals  are  held  together  by  second-order  (van  der  Waals) 
interactions,  it  is  clear  that  a  non- negligible  threebody  component  in  se¬ 
cond  order  of  perturbation  theory  may  play  a  decisive  role  for  the  stability 
of  such  solids. 

We  will  again  make  use  of  a  Gaussian  effective- electron  model,  i.  e. 
one  electron  per  atom  with  a  charge  distribution 


P  (r)  -  exp(-  p  ^r^)  ; 


(1) 


P  is  a  parameter  characteristic  for  the  rare  gas  atom  considered. 

Second- order  interactions.  Consider  a  triangle  (abc)  of  atoms,  and 
three  effective  electrons  1,  2,  3.  We  need  not  include  spin  variables  in 
the  wave  function,  so  that  we  write  in  zero  order  (single  Slater  determinant) 


.1/ 

*o  “  ^  ■  ^^abc)  ^  ^  ^  (2)  9c  •  ••  • 


where  the  overlap  integral  ^  is  given  by 


2  2  2  2 
A  «A,+A  +A.  -  2a  1.  a  At.  . 

abc  ab  ac  be  ab  ac  be  ' 
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in  terms  of  overlap  Integrals  between  pairs  of  atoms. 
For  the  second-order  energy  we  have  to  evaluate 


r*  ^  abcW^  abc^Ho  _ 

2  “  2_  — - 

K/to  ° 


abc 


<H' 


abc 


2 

>)  >....(3) 


where  is  an  "average  excitation  energy"  determined  by  the  averaging 
procedure,  is  the  perturbation  Hamiltonian  for  the  system  of  three 

Gaussian  atoms,  X  numbers  the  excited  states  of  the  system  (energy  ), 
and  where  is  the  unperturbed  ground  state  energy.  The  brackets  in¬ 
dicate  an  expectation  value  for  the  unperturbed  ground  state. 


Since  can  be  taken  directly  from  the  first-order  calculations 

the  only  unknown  quantity  is  .  We  denote  the  sum  of  second- 

order  energies  between  the  three  isolated  pairs  of  atoms  by 


E  ■  E  +  E  +  E 

*2  ‘^2  (ab)  ^  ^2  (ac)  *  ‘^2  (be)' 


and  form  the  difference  between  E^  and 


•  ®2  -  ®2‘°’  • 


This  difference  gives  the  threebody  component  of  Eg,  which  we  will  evalu¬ 
ate  in  relative  magnitude 


AEg/Egt®^  .-  (Eg  -  . 


2  “2 


(4) 
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The  expression  for  E^'  is  similar  to  that  for  E^,  with  replaced 

by  H'  ,  .  It  should  be  noted  that  the  average  excitation  energies  E.  occur- 
ring  in  the  expression  for  Eg  and  Eg  ’^ab)  necessarily  the  same 

quantity.  Their  difference  can,  however,  be  ignored  for  the  present  pur¬ 
poses  In  forming,  then,  the  ratio  A  Eg/ Eg  the  quantities  E^^  cancel 

to  a  sufficient  degree  of  accuracy. 

The  next  step  consists  in  substituting  the  explicit  expression  for  the 
perturbation  Hamilton  which,  of  course,  is  simply  the  sum  of  the 

three  pair- contributions,  i.  e. 


H'  +  H*  +  H*. 
ab  ac  be 


together  with  the  determinantal  form  for  into  eq.  (4),  and  evaluating 
the  various  expectation  values  which  occur.  It  is  essential  to  avoid  the  use 
of  multipole  expansions;  therefore  we  retain  In  its  exact  form. 

The  resulting  equation  for  Eg^®^  is  relatively  simple.  Let  us  intro¬ 
duce  the  notation 


V/K  <■»’' b' '2'  "'J  "’r  <*’2 

and  (5) 

'*•'1  ‘>^2 


C 


5  - 


for  direct  and  exchange  two-center  Integrals,  respectively,  with  similar 
expressions  for  and  ®bc* 


E*  E„ 
Av  2 


(o) 


D  ^  -  E  .  A  . 
ab  ab  ab 


(ab) 


1  -  A 


ab 


so  that 


-  E  .  E,^°) 
Av  2 


°ab  ■  ®ab  ^ab 


‘ab 


D  -  E  A 
ac  ac  ac 


1  -  A. 


ac 


°bc"  \c^c 


1  -  A 


be 


<  H' 


ac 


y  +  <  H' 


be 


(6) 
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On  the  other  hand,  the  formal  expression  for  <  >  is  already 

of  considerable  complexity.  The  square  of  contains  36  contributions, 
which  leads,  combined  with  different  integrals.  If  we  now 

decompose  into  squares  and  double  products. 


H'  K  ^  -  H'  ,  ^  +  H'  ^  +  H'  ^  +  2H'  .  H'  „  +  2H'  .  H'  ^ 
abc  ab  ac  be  ab  ac  ab  be 


t  2H'  H'  , 
ac  be* 


then  there  are  in  total  126  integrals  to  be  evaluated,  many  of  three- center 
type. 


Fortunately,  it  appears  that  all  these  contributions  to  the  second- 
order  energy,  and  also  the  expression  for  (eq.  6),  can  be  written 

as  linear  combinations  of  the  following  five  basic  integrals 

2  *  2 

A(x)  »  |e'*  J  e*  dt, 

o 


B(x) 


erf  .  X 

X 


C(x) 


2  , 
e'^  (erf.  -H—  y  du 

i2 


(7) 


F(x) 


I  I  e- 1  +  (u  +  ]  I  dt  j  du. 


where  x  is  a  side  of  the  triangle  (abc),  and 


00  2  71  71 

I  (u,  V,  e )  “  J  J  J 

000 


2ux  sin  a  sin 


2vxsina 


sin  (P  +8)] 


•  erf  [  (u^  +  -  2ux  sin  a  slnp  )  erf  [  (v  ^  +  x^  -  2vxsinasin  (p  +  0))  ^]* 


2  -x‘ 
X  e 


sinar  dad^  dx  , 


(7) 


where  u  and  v  are  two  sides  of  the  triangle  and  8  the  angle  between  them. 

In  principle,  the  basic  integrals  (7)  must  be  evaluated  by  electronic 
computation.  It  is  fortunate,  and  rather  surprising,  to  discover  that  asymp¬ 
totic  series  expansions  for  these  functions  can  be  given  with  excellent  pre- 
cision.  The  precision  is  of  the  order  of  10  for  small  values  of  the  variab¬ 
les;  this  corresponds  to  the  case  of  an  equilateral  triangle  of  the  heaviest 
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rare  gas  atoms.  The  precision  is  10  for  larger  values  of  the  variables 

•6 

(triangles  with  larger  opening,  and  lighter  atoms).  An  accuracy  of  lO’  is 
sufficient  for  the  present  problem. 

The  asymptotic  development  makes  it  possible  to  evaluate  the  final 
result  in  analytic  form;  the  expression  is  of  the  type  of  a  sum  of  exponen¬ 
tially  decreasing  functions  of  the  three  sides  of  the  triangle  considered, 
plus  a  sum  of  inverse  powers  of  these  distances. 

We  consider  here,  as  in  the  first-order  calculations,  triangles  form¬ 
ed  by  a  central  atom  and  two  of  its  twelve  nearest  neighbors;  0  measures 
the  angle  formed  by  the  two  sides  of  the  triangle  originating  at  the  central 
atom  (see  ref.  1  for  details).  Detailed  results  were  obtained  for  an  equi¬ 
lateral  triangle,  a  120°-conflguratlon  and  a  linear  symmetric  array  of 
argon  atoms;  approximate  values  were  already  determined  for  the  other 
important  triangles  0  =  110°  and  8  =  146°  of  the  hep- lattice.  In  addition, 
the  behavior  of  as  a  function  of  8  for  the  other  heavy  rare  gas 

atoms  was  established. 
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There  appear  to  be  two  Important  contributions  to  the  relative  three- 
body  component 

i)  effects  involving  diatomic  exchange  only,  i.  e.  termsarising 
from  exchange  of  electrons  between  two  of  the  three  atoms  on 
the  triangle; 

ii)  effects  due  to  triatomic  exchange,  i.  e.  those  involving  exchange 
of  electrons  between  all  three  atoms  on  the  triangle. 

In  the  following  Fig.  1  we  plot  these  two  components  as 

a  function  of  the  angle  0,  for  a  triplet  of  argon  atoms.  The  resulting  total 
effect  is  also  given.  (Fig.  1) 

It  is  seen  that  the  diatomic  exchange  effect  starts  from  a  rather 
high  value  at  60°  (0.  28),  decreases  rapidly  until  0  11  0°,  and  then  levels 

off  very  markedly  (0.  083  for  0  =  120°  and  0.  070  for  9  =  180°),  remaining 
positive.  Since  is  negative,  we  see  that  diatomic  exchange  favors 

the  formation  of  triangles  with  small  opening.  On  the  other  hand,  the  tri¬ 
atomic  exchange  effect  has  opposite  sign;  it  starts  at  a  very  high  negative 
value  (-0.  49)  at  0  =  60°,  decreases  much  more  rapidly  than  the  diatomic 

exchange  until  0  **  120°  (-0.  028)  and  then  remains  almost  constant  (-0.  01 

o  • 

at  0  =  180  ).  Consequently,  triatomic  exchange  favors  triangles  with  large 

opening. 

The  total  result  for  i.  e,  the  sum  of  diatomic  and  tri¬ 

atomic  exchanges. has  a  very  surprising  0-dependence  ;  the  curve  is  practi¬ 
cally  identical  with  the  one  obtained  in  the  first-order  calculations  (cf.  Fig.  1 
of  ref.  1).  The  resulting  value  is  approximately  -0.  20  at  60°,  increases 
rapidly  until  0  ^120°  and  remains  nearly  the  same  for  higher  values  of  0. 

At  0  =  110°  the  total  value  is  close  to  zero.  From  the  analytic  form  of  the 
final  result  we  conclude  further  that  AEg/E^^®^  for  neon,  krypton  and  xenon 
behave  in  the  same  way  as  for  the  first-order  effect. 
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Relative  second-order,  threebody  interaction 
for  triangles  a  =  b  =>  1,  c  (units  of  nearest  neighbor  distance) 
as  a  function  of  the  angle  0  between  the  sides  a  and  b.  The 
curves  I  and  II  denote  diatomic  and  triatomic  exchange  con¬ 
tributions,  respectively;  (I  +  II)  is  the  total  second-order 
effect. 
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We  have  now,  in  first-  and  second-orders  of  perturbation  theory, 
two  families  of  characteristic  curves  for  the  relative  threebody  inter¬ 
actions  between  atoms  of  the  heavy  rare  gases  which  are  very  close  indeed. 
Since  these  curves  are  all  of  the  same  type,  the  precise  values  of  the  di¬ 
mensionless  parameter  p  R  (R  is  the  nearest  neighbor  distance  in  the  lat¬ 
tice)  do  not  matter.  The  following  stability  considerations  apply  for  neon, 
argon,  krypton  and  xenon  as  well. 

Let  us  consider  the  resulting  second-order  curve  for  argon.  It 
was  shown  in  ref.  1  that  a  curve  of  this  type  favors  the  hcp-structure  in 
first  order.  Since,  however,  we  now  have  that  the  second-order  energy 
^2^°^  negative,  the  same  curve  stabilizes  the  fcc-latttce  in  second 
order.  Since  the  lattice  energy  is  negative,  the  second-order  effect  will 
be  more  important  than  that  of  first  order,  so  that  the  face-centered 
cubic  structure  is  the  more  stable  one  for  crystals  of  the  heavy  rare  gas 
atoms. 

To  calculate  this  net-effect  more  precisely,  we  determine  the  ratio 

AEj  + 


Consider  first  an  equilateral  triangle  of  (argon)  atoms,  placed  at  nearest- 
neighbor  distances  from  each  other.  In  this  configuration  we  have  ^2^°^  ^ 

-  2Ej^°\  Further,  since  we  found  that  the  first-  and  second-order  relative 
effects  are  very  nearly  the  same,  we  put  aEj/Ej^°^  ■  The 

above  ratio  becomes  now 


Ej^°4  e^°^ 


-  AE 
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E 

^1 


^^2  *  ^^2 


+  2 


-  10  - 


80  that  the  total  relative  threebody  effect  is  given  directly  by  the  second- 
order  contribution.  For  values  of  6  different  from  60°  the  ratio  -E2^°VEj^°^ 
for  the  triangles  considered  becomes  larger  than  2.  This  change  is,  how¬ 
ever,  very  slow  so  that  no  appreciable  error  is  involved  in  assuming  a 
value  2  for  all  these  configurations.  Further  we  find  that  AE^~-2 
the  same  relation  as  there  exists  between  E2^°^  and  Ej^°\ 

We  note  also  that  the  second-order  effect  agrees  in  sign  with  the 

4151 

Teller- Axilrod  triple- dipole  contribution  '  Some  of  the  general  charac¬ 
teristics  of  this  effect  can,  therefore,  be  taken  over  directly  from  the 
third-order  triple-dipole  interaction.  The  difference  between  the  hep-  and 
fcc-structures  is,  however,  at  least  an  order  of  magnitude  larger  than  the 
third-order  effect.  Further,  whereas  the  triple-dipole  effect  averages  out 
to  zero  at  distances  large  compared  to  that  between  nearest  neighbors 
because  of  its  angular  dependence,  the  present  threebody  interactions  are 
of  exchange  type  and,  therefore,  per  se  of  short  range. 

In  addition  to  explaining  the  stability  of  the  face- centered  cubic  con¬ 
figuration  for  solids  of  the  heavy  rare  gases  the  theory  predicts,  in  par¬ 
ticular,  that  the  interaction  between  atoms  in  an  equilateral  triangular 
array  will  be  20  percent  less  negative  than  the  sum  of  its  pair-potentials. 

In  this  connection,  a  very  recent  investigation  of  the  specific  heat  of  solid 
argon  (at  constant  volume)  by  Foreman  and  Lidiard^^  is  of  great  interest. 
For  their  calculation,  these  authors  used  the  enharmonic  Einstein  model; 
from  their  results  they  estimated  the  free  energy  of  formation  of  a  va¬ 
cancy.  There  is  a  discrepancy  between  these  results  and  experimental 
values^^"^\  Foreman  and  Lidiard  suggest  after  careful  analysis  that  the 
discrepancy  can  be  due  to  the  effect  of  threebody  interactions  and  estimate 
a  relative  effect  of  about  -25  %  for  an  equilateral  triangle  of  argon  atoms. 
Our  theoretical  result  is  in  agreement  with  their  estimate. 
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The  established  stability  of  the  cubic  crystal  structures  for  the 
heavy  rare  gas  atoms  is,  so  far,  only  valid  at  0°K.  At  non- zero  temper¬ 
atures  the  free  energy  determines  stability.  However,  Barron  and  Domb^^^ 
found  that  with  increasing  temperature  the  tendency  shifts  towards  the 
cubic  lattice  for  pair-potentials,  although  a  transition  from  hexagonal  to 
cubic  occurs  only  in  the  neighborhood  of  the  melting  points  of  the  solids. 
Since  the  threebody  component  is  of  short  "range,  with  a  relatively  large 
difference  between  the  two  lattices,  a  shift  towards  the  hexagonal  lattice 
at  such  low  temperature  is  unlikely. 

It  may  possibly  also  be  anticipated  that  the  unexplained  difference 
in  crystal  structures  occurring  among  the  alkalihalides  (for  example  NaCl 
and  KCl  on  one  hand,  CsCl  on  the  other  hand)  will  find  a  similar  explan¬ 
ation.  We  note  that  the  corresponding  ions  are  isoelectronic  with  the  rare 
gas  atoms.  Preliminary  investigations  are  in  progress  on  this  subject. 

We  emphasize  again  that  the  efficiency  of  the  Gaussian  model  in 
this  particular  case  may  be  mainly  due  to  the  fact  that  we  are  looking  for 
the  explanation  of  a  phenomenon  common  to  a  number  of  different  solids. 

As  a  consequence,  we  can  base  the  calculations  on  a  model  which  is  more 
primitive  than  that  required  for  a  specific  property  of  a  specific  atomic 
species. 

Detailed  results  on  all  different  triangular  configurations  will  be 
published  elsewhere.  Finally,  we  wish  to  thank  Mr  R.  Lenoir  and  his 
associates  of  the  Computing  Section  for  their  assistance  with  the  numeric¬ 
al  calculations. 


Genfeve,  le  30.1.1963 
Lj/SZ/lf  190-18 
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IMPLICATION  OF  RESULTS  FOR  POSSIBLE  FUTURE  WORK 


With  the  completion  of  the  analysis  concerning  the  stability  problem 
of  heavy  rare  gas  solids,  one  of  the  principal  aims  of  the  research  project 
has  been  achieved.  We  hope  that  these  results  may  be  classified  as  a  break¬ 
through  in  solid-state  theory,  since  they  concern  the  first  absolute  stabili¬ 
ty  proof  for  any  crystal  structure  occurring  in  nature. 

In  our  opinion,  we  are  also  approaching  the  final  solution  to  a  different 
major  problem  in  crystal  physics  :  the  explanation  for  the  crystal  structures 
of  the  heavy  alkali  halides  (e.  g. ,  CsCl  compared  with  NaCl  and  KCl)  and  the 
transitions  (induced  by  pressure  or  temperature)  occurring  between  these 
two  crystal  forms.  There  is  an  apparent  close  connection  between  the  alkali 
halides  and  the  heavy  rare  gas  atoms  :  the  ions  of  the  alkali  halides  are  iso- 
electronic  with  the  rare  gas  atoms  (Na  ,  F  with  Ne;  K  ,  Cl  with  A;  Rb  , 

Br  with  Kr;  Cs  ,  I  with  Xe).  It  has  been  found  in  the  literature  that  the 
Madelung  energy,  plus  pair- potentials  (van  der  Waals  interactions),  favors 
the  NaCl- configuration  by  as  much  as  one  percent  and  the  conclusion  has 
been  that  the  CsCl-structure  should  not  occur  at  all  (Cf.  M.  Born  and  K. 
Huang,  "Dynamical  Theory  of  Crystal  Lattices"  (Oxford,  (1954),  Chap.  III). 
Since  in  this  case  the  difference  in  pair-energy  is  so  large,  we  understand 
already  why  the  transition  from  one  structure  to  the  other  with  the  alkali 
halides  occurs  only  at  the  end  of  the  isoelectronic  series. 

Now  that  direct  experimental  evidence  for  manybody  interactions  in 
rare  gas  solids  also  appears  to  have  been  established  from  the  free  energy 
of  vacancy  formation,  we  can  envisage  that  the  elastic  constants,  for 
example,  of  these  solids  will  be  measured  in  the  near  future.  The  theo¬ 
retical  analysis  of  the  deviations  from  the  Cauchy-relations,  envisaged  in 
original  project,  falls  precisely  in  line  with  these  recent  developments. 
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It  is  interesting  to  observe  that  the  occurrence  of  manybody  inter¬ 
actions  in  heavy  rare  gas  solids  may  be  seen  as  a  direct  parallel  to  the 
recently  discovered  chemical  reactivity  of  xenon.  In  our  opinion,  very  sig¬ 
nificant  results  will  be  obtained  in  the  near  future  on  different  phases  of 
chemistry  and  solid-state  properties  of  closed- shell  atoms  and  molecules. 

The  interest  of  band  theoretical  calculations  for  such  solids,  in  con¬ 
junction  with  an  investigation  of  the  relation  between  symmetry  and  band 
structure  (and,  consequently,  electrical  and  optical  properties),  which 
we  have  already  initiated,  may  increase  considerably  as  a  result  of  these 
developments. 


PUBLICATIONS  AND  CONFERENCES 
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